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Abstract 



We briefly describe the simplest class of affine theories of gravity in multidimensional space- 
times with symmetric connections and their reductions to two-dimensional dilaton - vecton 
gravity field theories (DVG). The distinctive feature of these theories is the presence of an 
absolutely neutral massive (or tachyonic) vector field [vecton) with essentially nonlinear cou- 
pling to the dilaton gravity (DG). We show that in DVG the vecton field can be consistently 
replaced by an effectively massive scalar field [scalaron] with an unusual coupling to dilaton 
gravity. With this vecton - scalaron duality, one can use methods and results of the standard 
DG coupled to usual scalars (DGS) in more complex dilaton - scalaron gravity theories (DSG) 
equivalent to DVG. We present the DVG models derived by reductions of multidimensional 
affine theories and obtain one-dimensional dynamical systems simultaneously describing cos- 
mological and static states in any gauge. Our approach is fully applicable to studying static 
and cosmological solutions in multidimensional theories as well as in general one-dimensional 
\& . DGS models. We focus on global properties of the models, look for integrals and analyze the 

structure of the solution spaces. In integrable cases, it can be usefully visualized by drawing 
a 'topological portrait' resembling phase portraits of dynamical systems and simply exposing 
global properties of static and cosmological solutions, including horizons, singularities, etc. For 
analytic approximations we also propose an integral equation well suited for iterations. 

1 Introduction 



The present observational data strongly suggest that Einstein's gravity must be modified. The 
combination of data on dark energy (DE) and growing evidence for some sort of inflation have 
generated a wide spectrum of such modifications (see, e.g., [1] - [10]). Superstring ideas suggested 
natural modifications but, in view of the serious mathematical problems of the present string 
theory, strictly definite predictions about concrete modifications of gravity are not yet available. 
Moreover, the phenomenon of dark energy was not predicted by string theory and its origin in 
the stringy framework proved to be rather difficult to uncover and understand. The problem 
of dark energy in string theory looks very deep and is related to many other complex issues of 
(non-existing) quantum cosmology but it also led to some beautiful and exciting speculations, like 
eternal inflation and multiverse (see, e.g., [11] and [12]). On the other hand, if we modestly try 
first to find a natural place for DE in classical cosmological models, which inevitably are essentially 
nonlinear and non-integrable, we better return to recalling the origin of general relativity and look 
for some options abandoned or not found by its creators. 
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For this reason, simpler modifications of gravity that affect only the gravitational sector are 
also popular, e.g., [13] - [16]. In essence, these modifications reduce to the standard Einstein 
gravity supplemented by some number of scalar bosons (the first example of such a modification 
was the old Jordan - Brans - Dicke theory). The main problem of this approach is that the 
origin of these scalar bosons is quite obscure, and there is no theoretical principle regulating their 
coupling to gravity. Of course, there exist some phenomenological and theoretical constraints but 
in general the spectrum of these models is too wide. 1 The modification proposed and studied in 
[18] - [22] satisfy some principles of geometric nature based on Einstein's idea (1923) 2 to formulate 
the gravity theory in a non-Riemannian space with a symmetric connection determined by a 
special variational principle involving a 'geometric' Lagrangian. This Lagrangian is assumed to be 
a function of the generalized Ricci curvature tensor as well as of other fundamental tensors and 
is varied in the connection coefficients. A new interpretation and generalization of this approach 
was developed in [19] - [21] for an arbitrary space-time dimension D. 

The connection coefficients define symmetric and antisymmetric parts of the Ricci tensor (sij 
and ciij) and a new vector a«. Supposing that in the pure geometry there is no dimensional 
fundamental constants (except the velocity of light relating space to time) we choose geometric 
Lagrangians giving dimensionless geometric action. The geometric variational principle puts fur- 
ther bounds on the geometry and, in particular, relates Oj to a^. To define a metric tensor we 
have to introduce some dimensional constant. Then we can find a physical Lagrangian depend- 
ing on this dimensional constant and on some dimensionless parameters. Thus obtained theory 
supplements the standard general relativity with dark energy (the cosmological term, in the limit 
a* = a,ij = 0), neutral massive (or tachyonic) vector field proportional to cij (vecton) and, after 
dimensional reductions to D = 4, with (D — 4) massive (or tachyonic) scalar fields. 

Further properties of the theory depend on a concrete choice of the geometric Lagrangian. The 
resulting physical theory can be described by the corresponding effective Lagrangian that depends 
on the metric, vecton and scalar fields. In the geometric theory, there are no dimensional constants, 
while any fundamental tensor has dimension of some power of length (assuming c = 1). We 
proposed a class of geometric Lagrangian densities depending on the symmetric and antisymmetric 
parts of the generalized Ricci tensor, , , and on a fundamental vector obtained by contracting 
the connection coefficients , i.e., a% ~ Y^ k . By requiring the geometric 'action' (the space-time 
integral of the geometric Lagrangian density) to be dimensionless we can enumerate all possible 
actions, see [19]. The corresponding Lagrangian may be defined as the square root of an arbitrary 
linear combination of these densities or a linear combination of their square roots. 

The most natural density of this sort in any dimension is the square root of det(s;.,- + fa^-), 
where I is a number. 3 . The effective physical Lagrangian is the sum of the standard Einstein term, 
the vecton mass term, and the term proportional to det(gij + Ifij) to the power v = 1/(D — 2), 
where g±j and f\j are the metric and the vecton field tensors conjugate to and ctij , 4 I is a 
parameter related to I. The last term has the dimensional multiplier, which in the limit of small 
field fij produces the cosmological constant. For Z) = 4 we therefore have the term first introduced 
by Einstein but now usually called the Born - Infeld or brane Lagrangian. For D = 3 we have the 
Einstein - Proca theory, which is very interesting for studies of nontrivial space topologies. 

The simplest dimensional reductions from D > 4 to D = 4 produce (D — 4) scalar fields and 
thus the complete theory is rather complex, even at the classical level. Its spherically symmetric 
sector is described by a much simpler (l+l)-dimensional dilaton gravity coupled to one massive 
vector and to several scalar fields. This dilaton gravity coupled to the vecton and massive scalars 

1 Restricting consideration to homogenous cosmologies, one can find that in dimensionally reduced supergravity 
(superstring) theory there emerge many massless scalar bosons that couple to gravity only, see, e.g., [17], [2]. 
2 References to his papers as well as to related work of Weyl and Eddington can be found in [18] - [21]. 
3 Einstein used as the Lagrangian Eddington's scalar density 1 det r\j \ , where r»j = s;j + aij . 
4 This unusual construction introduced by A. Einstein is described and generalized in [18] - [21], see also Appendix. 



2 



as well as its further reductions to one-dimensional 'cosmological' and 'static' theories were first 
formulated in [18] - [21], and here we begin systematic studies of their general properties and 
solutions. These studies can be somewhat simplified by transforming the vecton field into a new 
massive scalar field, which is possible (on the mass shell) in the two-dimensional reductions. 5 As 
distinct from the normal scalar matter fields, the scalaron has a different coupling to gravity and 
may have abnormal signs of the kinetic term (phantom) or mass term (tachyon). Nevertheless, 
some general methods used in DGS models can be successfully applied to DG coupled to abnormal 
scalar fields. Thus, the main purposes of the present study are: to essentially modernize available 
tools for constructing various approximation schemes, to understand most important qualitative 
features of the global solutions, and to look for possible integrable approximation in simple one- 
dimensional reductions of the general theory. 

In this paper, we consider the simplest geometric Lagrangian, 



-Cgcom = y-det(sy +/<%•) = \/-A s , (1) 

where the minus sign is taken because det(sy) < (due to the local Lorentz invariance) and we 
naturally assume that the same is true for det(sjj + Tciij) (to reproduce Einstein's general relativity 
with the cosmological constant in the limit I — > 0). Following the steps of Ref. [19] or using results 
described in Appendix 7.1 one can derive the corresponding physical Lagrangian 



£phys — \/ ~9 



-2A [det(5j + lfi)\ u + R(g) - m 2 g ij a iaj 



„=!/(£> -2), 



which should be varied with respect to the metric and the vector field; m? is a parameter depending 
on the chosen model for affine geometry and on D (see [18] - [21]). This parameter can be positive 
or negative and we often use notation m 2 = fi. When the vecton field vanishes, we have the 
standard Einstein gravity with the cosmological constant. For dimensional reductions from D > 5 
to D = 4, we can obtain the Lagrangian describing the vecton a«, ~ diaj—djai and (D—4) scalar 
fields a/,, k = 4, .., D. Note that, for D = 3, the Lagrangian (2) is bilinear in the vecton field, and 
in the approximation ai = it gives the three-dimensional gravity with the cosmological term. 
The three-dimensional gravity was studied by many authors (see, e.g., [23], [24] and references 
therein) and this may significantly simplify the study of solutions of the new theory. 

We first consider the simplest spherical dimensional reduction from the .D-dimensional the- 
ory (2) to the two-dimensional dilaton-vecton gravity (DVG) and then further reduce it to one- 
dimensional static or cosmological theories. The next step is considering cylindrical reductions of 
the D-dimensional theory to more complex two-dimensional dilaton gravity theories coupled to 
several scalar fields. For simplicity, we consider this dimensional reduction only in the dimensions 
four and three. As is well known, in the D = 4 case, in addition to the dilaton, there appear 
'geometric' cr-model scalar fields and vector pure gauge fields that produce an effective poten- 
tial introduced in [25]. These two-dimensional scalars look like scalar 'matter' fields obtained by 
dimensional reductions from higher dimensional supergravity theory. Like those fields, they are 
classically massless. However, in general, the cr-model coupling of the scalar fields is supplemented 
by the mentioned effective potential depending on the scalar fields and on the 'charges' of the 'ge- 
ometric' pure gauge fields. This theory is rather complex and not integrable even after reduction 
to dimensions 1+0 or 0+1 (see [25]). Adding the vecton coupling makes the theory even more 
complex and in this paper we consider only the simplest variants. 



5 In addition to some formal simplifications, the transformation may shed a new light on some theoretical problems, 
especially, in cosmology. Indeed, new cosmological theories operate with several exotic scalar fields (inflations, 
phantoms, tachyons), which are usually introduced ad hoc. In the models considered here, such 'particles' appear 
as effective fields in a consistent theoretical framework. 



3 



Below, we mainly discuss low-dimensional theories (D = 1, 2, 3, 4) that allow us to consistently 
treat some solutions of realistic higher dimensional theories. We mostly concentrate on mathe- 
matical problems and do not pay much attention to physics meaning of the obtained solutions. 
Note also that in the context of modern ideas on inflation, multiverse, etc., the main parameters 
of a fundamental theory of gravity cannot be theoretically determined. In particular, we do not 
know the sign and the magnitude of the cosmological constant and of other parameters giving a 
scale of length. Even the dimension and signature of the space-time should be considered as a free 
parameters that can be estimated in the context of a concrete scenario of the multiverse evolution 
or by anthropic considerations. In practice, this means that we should regard theories with any 
parameters in any space-time dimension as equally interesting, at least, theoretically. 



2 Dimensional reductions of the generalized gravity theory 

Let us outline the main reductions of the model (2) in the dimensions D = 3, 4. 6 Due to natural 
space and time restrictions we only give an overview of main points. First, consider a rather 
general Lagrangian in the D-dimensional spherically symmetric case (x° = t, x 1 = r): 



ds l D = ds\ + ds l D _ 2 = gij dx 1 dx 3 + (p 2l/ d£L l D _ 2 , 



(3) 



where v = (D — 2) 1 . The standard spherical reduction of (2) gives the effective Lagrangian 7 the 
first three terms of which describe the standard spherically reduced Einstein gravity: 



'"D ~ V 9 



<pR{g) + k v <p'- 2v + W&) (Vv?) 2 + X(<p, P) - mW 



(4) 



Here ctj(i, r) has only two non- vanishing components do, a±, fij has just one independent component 
foi = a o,i — ai,o; the other notations are: a 2 = a^a* = g % 3 a^aj , f 2 = fijf tJ , k u = k(D — 2){D — 3), 
W(ip) ='(1 - v)/ip and, finally, 8 



X(tp,f 2 ) = -2Kip 



1 + ^A 2 f 2 



(5) 



where, the parameter A is related to dimensionless parameter I in (2) but [A] = L . 

Sometimes, it is convenient to transform away the dilaton kinetic term by using the Weyl 
transformation, which in our case is the following (w'(tp)/w(<p) = W{(p)): 



9ij=gijW (ip), w(Lp) = (p 



l-u 



S.2 



w z f , 



WSL 2 



(6) 



Applying this transformation to (4) and omitting the hats we find the transformed Lagrangian 



5(2) 
-D 



-9 



<pR(g) + K ~ 2Av" ( 1 + ^AV^-^f 2 



2 2 

m fa 



(7) 



When D = 3 we have v = 1, k v = 0, Weyl's transformation is trivial and the Lagrangian is 



,(2) 



'-9<P 



R(g) -2A- A*Af* -m 2 a 2 



(8) 



6 A careful consideration of the general successive dimensional reductions of Einstein's gravity from D — 4 to 
D — 1 is given in [26], [27]. For a review of the standard spherical and cylindrical reductions and solutions see [28]. 

7 We suppose that the fields <p, o< are dimensionless while [t] = [r] = L and thus [fij] = L _1 , [R] = [k v ] = [X] = 
[m 2 ] = L -2 . For more details on our dimensions see Appendix. 

8 The expression for the determinant in Eq.(2) in terms of fij for D = 4 written in [19] contains also the term of 
the fourth order in the vecton field. It is not difficult to see that in both spherical and cylindrical reductions (see 
the end of this Section) this term vanishes and thus Eq.(5) is valid. In fact, this is true in any dimension D. 
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These two-dimensional reductions are essentially simpler than their parent higher dimensional 
theories. In particular, we show that the massive vecton field theory can be transformed into 
a dilaton - scalaron gravity (DSG) model which is easier to analyze. Unfortunately, these DSG 
models and their further reductions to dimension one (static and cosmological reductions) are also 
essentially non-integrable. It is well known that the massless case, being a pure dilaton gravity, is 
classically integrable (see, e.g., [29] and reference therein). Having this in mind, we will attempt 
to find some additional integrals of motion in similar vecton models. 

The next simplified theory is obtained in a cylindrically symmetric case. We consider here 
only D = 3 and D = 4. The general cylindrical reduction was discussed in detail in [25] and here 
we only summarize the main results. The most general cylindrical Lagrangian can be derived by 
applying the general Kaluza reduction to D = 4. The corresponding metric may be written as 

ds\ = (gij + ip a mn <p™<p]) dx l dx j + 2p im dx^y 171 + ip a mn dy m dy n , (9) 

where i,j = 0,1, m,n = 2,3, all the metric coefficients depend only on the x-coordinates (t,r), 
and y m = (<p,z) are coordinates on the two-dimensional cylinder (torus). Note that p plays the 
role of the dilaton and o mn (det a mn = 1) is the so-called cr-field. The reduction of the Einstein 
part of the four-dimensional Lagrangian, ^J—g^ i?4, can be written as: 



J-4c — V 9 



2 



raj 



(10) 



where p^- = dip™ — djpf . The Abelian gauge fields p™ are not propagating and their contribution 
is usually neglected. In [25] we proposed to take them into account by solving their equations of 
motion and writing the corresponding effective potential (similarly to what we are doing below in 
the spherically symmetric vecton gravity). Introducing a convenient parametrization for a mn , 

CT 2 2 = e v cosh £, 0-33 = e~ v cosh £, cr 23 = cr 32 = sinh £ , 

one can exclude the gauge fields ip™ and derive the effective action 

(11) 



*-eff — V 9 



pR{g) + ^(V(^) 2 - |[(V£) 2 + (cosh£) 2 (Vr,) 2 ] + V eS (p,^, V ) 



Here the first three terms are standard, while the effective geometric potential, 

cosh£ 



2p 2 



Qle'" 1 -2Q l Q 2 t<m\ii + Q 2 2 e ri 



(12) 



depends on two new arbitrary real constants Q m , which may be called 'charges' of the Abelian 
geometric gauge fields p™ . This representation of the action is more convenient for writing the 
equations of motion, for further reductions to dimensions (1 + 0), and (0+1) as well as for analyzing 
special cases, such as Q1Q2 = 0, ^77 = 0. 9 

The static solutions of the theory (12) with Q1Q2 7^ have horizons while the exact solutions 
discussed in [25] for Q\ = Q2 = and non-vanishing u-fields £, 77 have no horizons at all, in 
accordance with the general theorem of papers [29]. An interesting special case can be obtained if 
we choose Q± ^ 0, Q2 = 0, £ = 0. Then, Lagrangian (11) gives a standard dilaton gravity coupled 
to the scalar field rj, with the potential depending both on the scalar field and dilaton. If Q\ 7^ 0, 
there exists a static solution with a horizon, which disappears when Q\ vanishes. Of course, the 
horizon also exists in pure dilaton gravity, when £ = 77 = 0. 



9 It is also closer to the original Einstein - Rosen equations for nonlinear gravitational waves, which can be 
obtained by putting Q\ — Q2 = and £ = 0. When Q1Q2 =fc 0, £ and r\ cannot be identically zero. 
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In [29] we studied in some detail the models with the potentials independent of the scalar. In 
this paper, we demonstrate that interesting results can be derived in more general models with 
'separable' potentials V(tp,ip) = vi(ip)v2(ip)- In particular, we show that one of the integrals of 
motion in the dilaton gravity coupled to massless scalars derived in [29] may exist also in some 
models with separable potentials that are of interest in the context of the present study. 

We see that the general cylindrical action is a very complex two-dimensional theory and even 
its one-dimensional reductions are rather complex and in general not integrable. Adding the vecton 
sector does not make it simpler and more tractable. Much more tractable is the three-dimensional 
cylindrical space-time. The metric can be obtained by the obvious reduction of (9), 

ds 2 = (gij + <Pi<Pj) dx t dx^ + 2</?j dx l dy + tp dy 2 , (13) 

and the corresponding Einstein Lagrangian is simply 

-4? = V-9 <p{R(g) - \ <Pij <p lJ } , <Pij = wj - w.i ■ (14) 

Using the equation of motion for 92^ and introducing the corresponding effective potential (see 
[25] and more general derivation below) we find the following two-dimensional dilaton gravity 

4ff = v 7 ^ {<t>R{g) - 8Q 2 <T 3 } <f> = • (15) 

As distinct from the cylindrical reduction of the four-dimensional pure Einstein theory correspond- 
ing to V e s = in (11), this theory has a horizon and can easily be integrated. 

It is not difficult to include the vecton part into the cylindrically symmetric Lagrangians. In 
fact, the terms —m 2 (pa. 2 , X(ip,f 2 ) in (5) are invariant and have the same form in any dimension. 
Therefore we can simply add the expressions (11), (14) to the gravitational Lagrangians also in 
the cylindrical case. Note however that there exist different cylindrical reductions of the vecton 
potential ctj. For example, unlike the spherical case, these fields may be nonzero for i = 0,1,2 
and correspondingly e 1 = f 01 = <9 ai - dia / 0, e 2 = /02 = d a 2 / 0, h 3 = f 12 = <9ia 2 / 0. 10 
We see that the component a 2 = a v of the vector field en behaves like an additional scalar field 
and, in addition to the two-dimensional vector field (ao,a±), we have up to three scalar matter 
fields, ci2,£,r], with a rather complex interaction to the dilaton gravity. This means that analyzing 
cylindrical solutions is more difficult than that of the spherical ones. Our consideration below is 
applicable to the simplest case of two potentials oq, a\ with one field /01 • 



3 Vector - scalar duality in two-dimensional vecton theory 

In the dimension D = 2 all fields (vector, spinor, ...) are practically equivalent to scalar ones. Such 
equivalence is widely known for massless Abelian gauge fields (see, e.g., [29] - [32] and references 
therein) . The aim of this Section is to establish a standard map of massive Abelian vector fields to 
scalar fields. We do not attempt to present the most general results in this direction and therefore 
restrict our consideration to a DVG coupled to scalar 'matter' fields. This includes all the models of 
the previous sections. Suppose that in place of the standard Abelian gauge field term, X(ip, iji) f 2 , 
the Lagrangian contains a more general coupling of the gauge field = di(ij — djdi to the dilaton 
and scalar fields, X(ip, ip; f 2 ), where f 2 = f % K Using the Weyl transformation we write a fairly 
general two-dimensional Lagrangian in the form 



cpR + V(<p, 1>) + X(<p, iP; f 2 ) + Z a (<p) a 2 + J2 Z(<P, V0(W0 5 



(16) 



Here e, = a oi , hi = Eijkdjk- In the diagonal metric gij = <?;<5;j the fourth-order term in the determinant from 
Eq.(2) is proportional to (eihi) 2 expj^ 2gk and is seen to vanish. This argument obviously works in any dimension. 
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where the kinetic dilaton term ~ ip' 1 (V(p) 2 is now absent while V and Z terms represent potentials 
and kinetic terms of scalar fields (like those in (11), (12)). We need not specify the number of 
the scalar matter fields and therefore omit the summation indices for the fields tj) and their Z- 
functions. For the fields, having positive kinetic energy, all the Z-functions in (16) are negative 
and usually proportional to ip ; we call them 'normal matter' fields. 

In [30], we considered the massless case, Z a = 0, and proposed to use instead of (16) the 
effective Lagrangian not containing the Abelian gauge fields (for a detailed proof see [31], [32]). In 
this paper, we extend this approach to the massive vecton and show that the vector field can be 
effectively replaced by a scalar one (scalaron) denoted by q(t,r). This means that the equations 
of motion defined by Lagrangian (16) can be obtained from the Lagrangian, 



,(2) 
-eff 



'~9 



ipR + vfr, vo + x cS (<p, v, g) + z q (v){v q ? + J2 YO(vy) 2 



(17) 



by applying a formally invertible transformation from the vecton fields ctj , to the scalaron fields 
diq, q. This can be constructed in a complete analogy to the massless case. 

First, consider the equations for the vecton fields derived from Lagrangian (16): 



dX 



di 



2 ' 



(18) 



At this point, it is convenient to rewrite the Lagrangian and equations in the LC coordinates, 

ds 2 = -4h(u,v)dudv , yf^g = 2h , fuv = a UtV -a VtU , -2i 2 = (f uv /h) 2 . (19) 

Before we pass to the LC Lagrangian we should derive the energy and momentum constraints 
using the general coordinates. Leaving only one extra field ip we have (see, e.g. [18], [29] - [34]): 



hdiidw/h) = Z a (ip)a 2 + Z(ip,ii)(drt) 2 , 



u, v. 



(20) 



The other equations of motion can be derived from the LC transformed Lagrangian, 



\cW = y dudv i n |^| + h [ V ^ ^ + X (<p, 4>, f 2 )] - Z a {<p) a u a v - Z(<p, d u ^ 3^ , (21) 



simply by variations in h, tp, f 2 , ip. The most important equations are the following: 



d u d v (p + h 



, dX df 2 



= d u d v (p + h 



X -2f 



2 dx 



df 2 



o, 



X(ip,^f 2 ) 

which allows us to define the effective potential X e g(ip,ip,q), and the vecton equations, 
d u {hr l X'f uv ) = -Z a (ip) a u , d v (h~ l X' f uv ) = Z a (ip) a v , 



(22) 



(23) 



that are equivalent to (18). Now we transform them into a standard scalar equation for a new 
scalar field q(u,v), which can be defined both in the LC and in general coordinates: 



2^r'xw^, i,j = o,i. 



(24) 



It is not difficult to understand that so defined scalar q is also invariant under the Weyl transfor- 
mation (6). With this definition, we may regard f 2 as a function of q, which can be derived by 
solving (24). Denoting the solution by z(q) and expressing f uv /h in terms of f 2 with the aid of 
(19) we find an important relation: 



2z(q) = -(q/X') 2 , 



X' = 



d 

d z(q) 



X[tp,ip;z(q)} . 



(25) 
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Returning to (22), we see that the expression in the square brackets, with f 2 replaced by z(q), 
is in fact our effective potential X e g(tp, ip, q). There exist several expression for it that allow one to 
derive the same equations of motion for the fields h, <p, ip, q. In [31], [32] we obtained two equivalent 
expressions for X e g by using (22) and (25)): 

X eS (cp, ^ q) = X[<p, V>; z(q)} - 2z(q)X' = X[<p, ^ z(q)} + q 2 /X'[<p, ^ z(q)} . (26) 

From (25) and (26) we also derived the most compact and beautiful representation: 

X eS {<p,i;,q) =X(ip,iP;z) + qy/=2t, (27) 

where we intentionally write z instead of z(q). The reason is that 

d q X cS (^,ij,q) = {X' -q/^Tz )d q z{q) + ^2z = yf=2z , (28) 

as follows from Eq.(25). Now it is easy to write the equation for q. Using (23), (24) we find: 

a u (u,v) = -Z~^(y)d u q(u,v) , a v (u,v) = Z' 1 ^) d v q(u,v) , (29) 

d^Z^duq) + OuiZ^dvq) = -f uv (q) = -hq/X' = -hd q X eS (ip,^,q) (30) 

The last equation for q(u, v) can be derived from Lagrangian (17) by taking Z q = Z' 1 . Usually, 
Z a = —m 2 tp and then the scalaron kinetic term has the unusual form, 

Z q (Vg) 2 = Z" 1 (Vqf = (Vqf/i-m^) . (31) 

It is not difficult to derive all the equation of motion using effective Lagrangian (17) with X e g and 
Z q given by (27), (31), and we list them here for further reference. In addition to the relations 
obtain above, one has to use the following easily checked identities: 

d h X eS (p,tp,q) = 0, d^Xcs = X(ip, tp, q) , d^X eS = X(cp,ip,q) . 

The first identity immediately follows from our definitions. In fact, it is a characteristic property 
of X e Q and thus can serve to define it. Two other identities are derived like Eq.(28). 



Let us define the total effective potential, 

U(v,tp,q) = V(<p,ij)+X cS &,^,q), (32) 
and first write the LC version of the effective Lagrangian and of the constraints 

1 (2) 

-C\4 = ipd u d v ln\h\ + hU(tp,ip,q) - Z q (ip) d u qd v q - Z(tp,ip) d u tp d v ip , (33) 

hdi(di<p/h) - Z q (cp) (d t qf - Z(^iP) (M) 2 = 0. i = u,v. (34) 
The other equations are the same as used in the dilaton gravity coupled to scalars: 11 

d u d v tp + h U(ip, ip, q) = 0, (35) 

d u (Z q d v q) + d v {Z q d u q) + hd q U{ip,^,q) = 0, (36) 

d u (z d v ip) + d v (z d u ip) = d^[-hu + z(ip,7p)d u ipd V 7p], (37) 

d u d v \n\h\ = d^[-hU + Z q (ip)d u qd v q + Z(ip,7p)d u ipd v iP]. (38) 



Below, we also use notation In \h\ = F, h = eexpF, where e = +1 for cosmological and —1 for static solutions. 
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One of the last three equations is satisfied by solutions of the remaining five. The energy- 
constraint, which is the sum the two constraints, plays a special role in the ADM Hamiltonian 
formulation [35]. We do not discuss here these subtleties that will be clarified in considering the 
one-dimensional reductions of this two-dimensional quasi-linear constrained system of the second 
order. The structure of this system allows one to find some integrable classes using simplifying 
assumptions about the potentials U and Z. For some multi-exponential potentials U and constant 
potentials Z, there exist integrable systems related to Liouville and Toda-Liouville ones (see [36], 
[37], [33]). The pure Liouville case was completely solved but, for the Toda-Liouville, it is difficult 
to find an exact analytical solution of the two-dimensional constraints, even in the simplest u\®su2 
case. This problem is easily solved in the one-dimensional (static or cosmological) case. 

For the spherical D = 3 and D = 4 Lagrangians in the Weyl frame (see (8) and (7), resp.), it 
is easy to derive the effective potential U = U/w(ip) and thus to find potential (32): 

U(<p,q) = U(<p, q) = -2Av9[l + g 2 /4A 2 AV 2 ], D = 3, (39) 

U{f,q) = v^U(cp,q) = -2A93[l + g 2 /A 2 AV ]^ + 2fc , D = 4. (40) 

We remind that the potential Z q is D-independent, Z q = — l/m 2 (p, see Eq.(31). 

As was argued in [18], even the one-dimensional cosmological reduction of the pure scalaron 
theory with the potential (39) is not integrable. Thus we concentrate on searching for approximate 
potentials that allow us to find a wide enough class of analytic solutions to reconstruct exact ones by 
iterations. For instance, potential (40) for large q is asymptotically very simple, in the standard 
frame we have U = —2q/\ + .... As will be shown shortly, this asymptotic formula is valid in 
dimensions D > 3. On the other hand, the behavior of the ^-dependent part of U for small q is 
given by (39), if we multiply the q 2 term by (D — 2). This universality allows us in what follows 
to regard (39), (40) the generic scalaron potentials. 



4 Cosmological and static reductions of scalaron theory 

In this paper, we consider only the simplest reductions of the two-dimensional theory ignoring 
one-dimensional waves studied in our previous work [30]- [33], [37]. The simplest reduced gauge 
fixed equations can be directly derived by supposing that in the LC equations (33) - (38) the fields 
h,ip,q,tp depend on one variable, which we denote r = (u + v). For the cosmological solutions 
this variable is identified with the time variable, r = t, while for the static states, including black 
holes (BH), it is the space variable, r = r. The only difference between the cosmological and static 
solutions is in the sign of the metric, h c > 0, h s < 0. 

In our study of black holes and cosmologies we use the general diagonal metric, 



ds 2 . 



3 2a(t,r) dr 2 _ ^r) df 2 



(41) 



Then, the static and cosmological reductions of our two-dimensional vecton theory (17) can be 
presented by the Lagrangian (with the standard Z(ip,ip) = —if) 



ip 2 - 2d e - - (1 - v) 



V + X(f* 



(42) 



Here we omit a possible dependence of V and X on tp and ip, denote Z a = —flip = m 2 , and e = ± . 
All the fields depend on r = t (e = +) or on r = r (e = — ). Finally, 



a + = ai (r) , a_ = ao (t) , a+ = a , a. 



d 



d 



7, OL e = -^0! e , a = — a. 
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We see that the cosmological and static Lagrangians essentially coincide, the only difference 
being in the sign of the potential term and of the metric exponents as well. If e = ±, the multiplier 
of the kinetic term, l e = e a ^ , is obviously a Lagrange multiplier varying of which produces the 
constraint equation, which is equivalent to vanishing of the Hamiltonian. In view of the implicit 
dependence of f 2 on l e , it is much more convenient to first employ the canonical formulation 
and then identify a proper Lagrange multiplier. Alternatively, one can follow the steps made in 
previous Section but the canonical approach is simpler both technically and conceptually. 

Let us construct the effective Hamiltonian and Lagrangian of the vecton sub-Lagrangian, 

C a = hX(f 2 ) - fop a 2 ; fi = e th ~ a) fi, h = ee a+ ~*. (43) 
Remembering our previous results we define the variable y and derive the momentum p a , 
1 , 2 „ 2 l2 /«\ 2 2 1 , dC a BX 

where the last equation implicitly defines y(p a ) ■ This allows us to write the partial Hamiltonian, 

ox- 



T-L a {Pa ,a) = p a a- £ a = -h 



+ flip a 2 = — hX e g + flip a 2 . (45) 



Here X e g essentially coincides with the effective potential derived in the previous section but now 
depends on ip,ip and y(p a )- Applying the simplest possible canonical transformation, 

p a ^-2q, a^p/2, X cS (p a ) => X cS (-2q) , (46) 

we see that the effective potentials (and the partial Lagrangians) obtained here and in Section 3 
are identical. Indeed, applying the canonical transformation to Hamiltonian (45), we find the 
partial Hamiltonian of the scalaron, % s {P-> q), and the corresponding Lagrangian, 

C,(p,q) = pq-H s (p,q) = hX eS + e< a ~^ ±- , (47) 

flip 

which can now be inserted into the complete Lagrangian instead of C a . 

Thus we have the complete DSG Lagrangian, now denoting l e = exp(a_ e — a+ e ), 



;2 r ;2 



^i? - 2a e ip - (1 - v) — + At- 

ip m^ip 



+ l e ee Ia *U{^,q), (48) 



where we also use notation (32) and n = m 2 . This form is more natural than (42). First, the 
dependence on the Lagrangian multiplier l e is simple and standard, the kinetic part is quadratic 
in generalized velocities and can be made diagonal by a redefinition of the Lagrangian multiplier 
and velocities. In addition, we are free to make a convenient gauge choice and to choose the Weyl 
frame. For example, by making the shift a e => a e — (1 — v) In ^/<p and redefining the potential by 
U =>■ ip u ~ 1 U we remove the third term in (48) and obtain the Lagrangian in Weyl's frame. 12 Then 
we can redefine l € ip = l e , introduce the new field £ = ip 2 and finally rewrite (48) in a simpler form, 

eC^ =l- l [^ 2 + m~ 2 q 2 -ia e ] + l e e e 2a < ^ /2 ~ l -U( VI , q, 1>) ■ (49) 



12 After this transformation, the theory defined by constrained Lagrangian (48) can be compared to the one- 
dimensional reduction of the theory defined by equations (33) - (38): we first derive the constraint, choose the 
LC metric by putting i e = f, and then derive all the equations of motion. Note that the general one-dimensional 
Lagrangian (48) gives more freedom in simplifying equations. For example, if there are only standard fields ip, we 
may choose l e = l e tp, £ = ]nip and then there will be no need to use Weyl's frame for simplifying the kinetic terms. 
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Before applying it to studies of cosmologies and horizons in the scalaron theory we discuss the 
effective scalaron potential, corresponding to the X-potential (5) in more detail. 

Using notation (44) it is not difficult to find that y{p a ) can be derived from the equation 

y = x{l-y 2 ) 1 -\ x = q/(-2v\Av), (50) 

where we choose the sign of q so that x > and therefore y > 0.. This is a polynomial equation 
of order (D — 3) for y 2 if D is even, and of order 2(D — 3) for y if D is odd. This observation is 
not very useful if D > 4, but some algebraic expressions for y(x) can be written if D = 5, 6. They 
are rather complex and difficult to use. However, it is not difficult to find good approximations 
for the solutions for small and large values of x: 

y = x[l- (l-u)x 2 + ...]; 1-y 2 = \x\~ a - (a/2)\x\~ 2a + ... , a = 1/(1 - v) . (51) 

Both expansions are applicable to all v > 1, the first is exact for v = 1 as y = x. 

To analyze the behavior of X e g{x) let us rewrite it in terms of the solution y(x) of Eq.(50) 



BX r 
X c z=X- y ^ = -2kv- 

oy y 



1 - (1 - 2v)y 2 



(52) 



This expression is convenient for deriving some exact solutions as well as for discussing their general 
properties for arbitrary D. If D = 3, then y = x, and we directly find the effective potential (39). 
For D = 4, x/y = a/1 + x 2 , and thus the potential coincides with (40). An interesting feature 
of this effective potential is its invariance under rotations in the plane (q,\Aip). This invariance 
exists only in D = 4. Unfortunately, the kinetic part does not respect this symmetry for any 
gauge choice. A different, simpler symmetry has the effective potential in D = 3, which is equal 
to (39) divided by <p = y/£ . Then Lagrangian (49) of the pure scalaron gravity, with h = e e 2cte 
and tp = 0, has a very simple form, 

eC^ =l- 1 [m- 2 q 2 -ia t ] - l e 2Ah [ 1 + q 2 /AX 2 A 2 £] , (53) 

in which it is easy to find the scale symmetry q >->■ C q, £ i— > C 2 h >->■ C 2 h . Such a symmetry 
itself does not allow to obtain and additional integral of motion. At most, it may signal of existence 
a special solution with a scaling-type dependence on q, £, h, but this is not our primary concern. 
Below, we will show that, in fact, the theory (53) has an additional integral, but first we should 
complete our discussion of the general properties of the theories (49). 

Applying asymptotic approximations (51) to (52) we prove the universality mentioned above: 



X cS = -2A<p [1 + q 2 /4v\ 2 A 2 <p 2 +0(x 4 )], X cS = -2^/ q y\ 2 + 0(\x\~ CT ). (54) 

The small q approximations is exact for v = 1, the large q approximation is independent both of 
if and v > 1, and the general structure of X c q is defined by a positive regular function, v u (x), 

X e s/(—2Aip) =v u (x); v v (0) = l; v u {x)~2vx, x — > oo . (55) 

Thus the effective potential in Lagrangian (49) is 

U e (<p,x) = r /2_1 ^(v^,(Z,^) = V U ~ 2 U= -2A^- 1 ^(x) + K<P~ {1+V \ (56) 

and the exact expressions for v = 1, 1/2 are given in (39), (40). 

Let us discuss properties of v u (x), which can be considered as the function of q 2 /(p 2 = q 2 /$, 
(it is expressed in terms of y(x) by Eq.(52)). Using (50) and (51) it is easy to find that 

d y y 1 ~ y2 >/n\ 1 U \ n 

^ = x l + (l-2^ >0 ' y( ° ) = 1 ' ^°°) = ' 
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v' v (x) = 2v y{x) ■ v' v {x) > , > . (57) 

Thus v u {x) is monotonic concave function having simple expansions for x <C 1, x S> 1: 

1 - i/ 



iv(x) = 1 + vx + 0{x ); v u {x) = 2vx 



1 + -7^ x ~ a + 0{x- 2a ) 



(58) 



With such a simple and regular potential U e , one might expect that at least qualitative behavior 
of the solutions of the theory (49) could be analyzed for small and large values of x. This is 
true if the theory is integrable. But it is not integrable, even in the simplest D = 3 case when 
vi(x) = 1 + x 2 , k\ = and thus U e = — 2A(1 + x 2 ) is linear in q 2 /£. The integrability of this and 
of similar but more general theories is discussed in next Section and in a separate paper. Here we 
consider the scalaron canonical equations and find one additional integral of motion. 

Consider the general pure scalaron theory in which we neglect the curvature term k u (p~( 1+u ^ : 

£W = J' 1 [ - £ a £ ] - l e 2Ah i/" 1 v u (x), n = m 2 . (59) 

We did not yet fix the gauge, T e is arbitrary. Varying in it we get the vanishing Hamiltonian, 13 

n[ 1] = l e H<U = l e [ np 2 q /A - Pa pt + 2Ah if"" 1 v v (x)]=0, (60) 

where we used the definition of the canonical momenta (here and below, a e is replaced by a), 

q=l e fip q /2, i = -T £ p a , a = -T e p£. (61) 

We may use the parametrization invariance of the theory to choose a most convenient evolution 
parameter r, for example, df = l e dr = £ dr = d^. 1A Such a parametrization is equivalent to gauge 
fixing plus a coordinate transformation. The standard LC gauge choice in our present notation 
is l e = 1, we mostly use the gauge l e <p = l e = 1 and keep the dot notation for d/dr. To further 
simplify the equations, consider the v = 1 case and denote 2hhv\ (x) = hu(q 2 /^). With this 
notation and with 4 = 1, the Hamiltonian system is given by equations (61) plus 

p q = -2hLo'{q 2 /£) q/i , Pi . = hJ {q 2 /i) q 2 /f , p a = -2hw . (62) 

It is easy to see that this Hamiltonian system has two integrals of motion. The first is the 
Hamiltonian (60). The second one follows from the obvious relation qp q + 2^p^ = 0: 

d/dr (qp q + 2^p ( ) = qp q + 2^p^ = -2hu = p a , (63) 

where we used (61), Hamiltonian constraint (60), and the last equation in (62). Integrating this 
relation gives the desired second integral. Thus we have two integrals: 

VPq ~ kPaPt + = ; qp q + 2£p i . - p a = c . (64) 

Knowledge of these integrals does not allow us to integrate our system but it is helpful in analyzing 
its properties. One may hope that, for some particular potentials u, a third integral can be found. 
For these reasons we call such systems partially integrable. The present observation suggests a 
more general approach to constructing system with additional integrals presented in next Section. 



13 Alternatively, we can use the standard Legendre transformation and write the Hamiltonian theory in the La- 
grangian form, £(q,q) — pq — lH(p,q), which is most convenient for systems with constraints (see, e.g., [38]). 
14 This choice is locally possible and is good as far as £ 7^ 0. Here we ignore such subtleties. 
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5 Integrals and integrability in simple cases 



Here we consider a general DGS with one scalar ip that may be a standard field or the scalaron: 



^-dgs — v y 



pR + Z(p)(V^) 2 + V(p,ij) 



(65) 



For the scalaron we have ip = q, Z = Z q = — l/(m<p) and the potentials are given above. For the 
standard scalar Z^ ~ —ip, but some of the results presented below are applicable to more general 
Z-functions. In our notation, negative signs of Z correspond to positive kinetic energies of the 
scalar fields but our classical consideration is fully applicable to both signs. The general model 
(65) with a general potential V is not integrable in any sense. One of the strong obstructions to 
integrability is the dependence of Z on ip, and the usual simplifying assumption is that the Z- 
functions are independent of ip. With this restriction, there exists a class of 'multi - exponential' 
potentials, for which the DGS theories with any number of scalar fields can be reduced to the 
Toda - Liouville systems and exactly solved. 15 For their static - cosmological reductions, analytic 
solutions were explicitly derived. Here we try to expand this class of the models. 

First, consider the case of the massless vecton. In the pure scalaron DG we have, in the 
massless limit, hdi(dipfh) = 0, from which we find h(u,v) = h(r), <p(u,v) = p{j), where r = 
a(u) + b(v). This means that the theory is automatically reduced to dimension one and we can 
use the Hamiltonian formulation of Section 4 (in the LC gauge). It is clear q is a constant, 
q = qo, as can also be seen from (23), (24). Now, using the above definitions of the momenta 
p a , , constraint p a P£ = hu(^), and equation p a = —2hco(^), all of which are valid for arbitrary 
potential u;(£), one can explicitly solve the equations of motion: 

p a = -2huj = -2p a p£ => 2ap a = p a => - £ = Pa = -c\h => ci 'h = -2hu , (66) 

where c\ is the integration constant. Then it follows the well known expression for h, 

c 2 1 — = 2uj(0, c 2 1 h = 2jd^u;(O = 2jdipV(ip)=N(ip)-N(ip ), (67) 

if we recall that the transformation of £ to <p induces the transformation of the potential, so that 
d£, w(0 = dipV(tp). Thus we find the explicit general solution of the general dilaton gravity, if we 
in addition solve the equation £ = c\h(^) defining the r-dependence. Now, in our special case of 
the dilaton - scalaron gravity corresponding to the massless vecton, the only trace of the vecton 
is the dependence of the potential on g , the 'charge' of the massless vector field, and on other 
parameters defining the theory. The general solution depends on the free independent parameters, 
(ci, ipo, To), 16 which can be hidden by rescaling of h and ip and by shifting r. Then the 'portrait' 
of the solution is a curve in the (h,ip) plane significantly depending on qo- The most important 
feature of this portrait is the structure of the set of the horizons. For us, the most interesting is 
the dependence of this portraits on the parameters of the parent DSG: qo,A,X,u. Characteristic 
properties of the portraits are: the number and structure of the horizons (simple, degenerate, 
singular) and the behavior of the solutions near horizons and singularities. A more interesting 
portrait of an integrable DGS is described in Appendix. 

For DGS (65), we can derive some convergent expansions of exact solutions near the horizons 
by applying the general approach proposed in [34]. Here we formulate its most compact form, 



15 This class includes all previously considered integrable two-dimensional DGS, which are reviewed in [36] (see 
also [39]). The first DGS of the Liouville type ('bi Liouville'), which generalizes the so-called Jackiw [40] and CGHS 
[41] models, was proposed and solved in papers [29], the results of which are essentially generalized here. 

16 The physical parameters are defined by the equation for horizons, N(ip) — N(tpo) = 0. When this equation has 
a unique non-degenerate solution, as in the case of the Schwarzschild black hole, N(ipo) can be related to its 'mass'; 
the case of many solutions, (p n , the physical interpretation is much reacher, see, e.g. [26], [42]. 
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which is convenient to use in general DGS with arbitrary Z, V {V^ / 0), and any number of 
scalars. For the sake of generality, we do not explicitly use the canonical formalism and thus 
somewhat change notation to the one resembling notation of [29]. We will also use the LC gauge 
and Weyl's frame having in mind that Eqs.(48)-(49) show how to return to the standard frame 
and general gauge. Thus the system of one-dimensional equations can be obtained by dimensional 
reduction of (34)-(38) if we suppose that all unknown functions depend on one variable r = u + v. 
Denoting F = Inh and introducing the new momentum-like variables \-> Vi Pi 

<p = X , Z(<p)iP = r ] , Z(<p)F=p, (68) 

we rewrite the main dynamical equations and the constraint as 

X + hV = 0, 2r ] + hV i , = 0, p + h(ZV) ip = 0, X p + hZV + t] 2 = , (69) 

where the lower indices <p, and ip denote the corresponding partial derivatives. It is well known that 
the second or the third equation can be omitted 17 and actually we have two independent second 
order differential equations and one constraint which is the first order equation. Therefore, to 
completely solve this system we must search for two additional constraints canonically commuting 
with the Hamiltonian and with each other. 18 

Now, recalling the previous section we change the variable r to £ that is defined by the relations: 
xdr = dip = Zd^. Then we rewrite equations (68)-(69) denoting the derivative d/d^ by the prime 
and introducing useful notation: U(^,tp) = Z(ip)V(ip, tp) , H(^) = h/x, and G(^) = rj/x- The 
main independent equations for x-> "*}■> ^ H now have a very compact form: 

i// = G, H' = -G 2 H, x '+UH = 0, 2r]'+U^H = 0. (70) 

The extended system contains two equations for p (see (68)- (69)): 

X F'- P = 0, p' + UsH = 0, X G' = UH(G-U^/2U), (71) 

where we add the explicit equation for G, which may be used instead of the last equation in (70). 

The system (70) is most convenient for deriving the solutions near horizons and in asymptotic 
regions and for studying their general properties. For example, a very important property of its 
solution is that (InH)' = —G 2 < 0. This property does not depend on the potential and is true 
for any number of scalar fields provided that their Z-functions are negative, as was first shown in 
[22]. Indeed, in this case the constraint equation can be written as: 

N 

<&'= (\nH)' = -Z J2Z-\0(rin/x)\ % = ^ = X ^/. (72) 

n=0 

For normal fields Z n < for all n, for some anomalous fields, like scalaron corresponding to 
tachyonic vecton, Z may be positive and then the sign of depends on the concrete solutions. In 
case of the same signs, Eq.(71) resembles the second law of thermodynamics and defines an 'arrow 
of time' for our dynamical system. If this is not true, the theorem is violated in a very specific 
way. It may be an interesting point for cosmological modeling. 

Our system of equations (70) has other interesting properties. The general solution of the first 
two equations can be written in terms of integrals of the function G(^): 

m = Vo + j^G, H(0 = H exp(-^G 2 ^j . (73) 

17 The full system (68)-(69) is evidently over-complete and may be called the extended system. It is useful when 
we look for integrals of motion. Its subsystems can be used for constructing some special solutions. 

18 As far as we are interested in the classical theory we usually will look for integrability in the Liouville sense. 



14 



Then, inserting these 'solutions' into the third and the forth equations and integrating them we 
can write one integral equation for G(£) instead of system (70): 

C<0 = * -(»" 5 £^)(*> -£<"*)"' • (74) 

where ip{£) and H{g) are given by Eq.(73). As was discussed in [29] and [34], the standard 
(regular and non-degenerate) horizon appears when Xo = Vo = 0. Then /i(£o) = while G(£o) = 
Ui/)(£o)/ 2U(£o) is finite if ^ 0. It follows that G, ip, H are finite and can be expanded in 
convergent series around £o if the potential £/"(£, if?) is analytic in a neighborhood of (£o> ^o)- 19 
When = there is the obvious integral of motion r\ = t/q. As can be seen from the above 
equations and was proved in [29], there is no horizon for r/o^O but, in the case of r/o = 0, we have 
G = and return to pure dilaton gravity with the solution (67) always having horizons. 

In simple cases, the integral equation can easily be reduced to a differential one. For example, if 
U = v(ip) and = 2g U, the integral equation can be reduced to the second-order differential 
equation, which is not integrable for arbitrary u(£) but is explicitly integrable if U$ = g\U. 
This result is quite natural as in this case there exist two additional integrals, r] = gx + t)q and 
p = gix + Po and therefore the most direct approach is to use the extended differential system. 
In non-integrable cases or when there is only one additional integral, the integral equation still 
can be a quite useful analytical tool, which we describe in a separate publication. Here, we only 
briefly outline a generalization of the approach of Ref.[29] to finding potentials U{^,ip) for which 
the extended differential system has additional integrals. 

Generalizing the approach of [29] and the above remarks about possible integrals of motion let 
us collect those equations which can generate such integrals: 

= p + UH + if/x = x' + UH = rj + Efy H/2 = p' + U t H = , (75) 

= = ntf-rf/x = Zp' + tU^H = ?p - p = 0, (76) 

where the first equation in (75) is the energy constraint, which we denote Eq and the next items 
in this chain of equations are denoted by , i = 1,...,7. Now, taking the sum J2o c i^i with 
C4 = C5 = cq = cj = Co we find that the solutions of equations (75) satisfy the identity 

[ciX+ c 2 r,+ c 3/ 9+ c 4 (?M+ £p)}' = -H[( Cl + c 4 )U+ c 2 U f /2+ c 3 U^+ 04(1^^/2 + ^U^)] . (77) 

Therefore, if the r.h.s. identically vanishes, the l.h.s generates the integral of motion, 

ciX + C 2V + c 3 p + c 4 (V>?7 + £p) = h . (78) 

This means that for the potentials U(£,ip), satisfying the partial differential equation 

(ci + c 4 ) U + c 2 Uip/2 + c 3 Ut+ c 4 (^ + f = , (79) 

there exist the corresponding integral of equations (75). All the above integrals can be obtained by 
applying this theorem. 20 The solution of Eq.(79) depends on an arbitrary function of one variable. 
Using this fact it is possible, in some simple cases, to derive one more integral. 

In the above example of the potential U = Uq exp(2gip + gi£) with two additional integrals we 
immediately find the equation for x-> 

X' = {9 2 + 9i) X + {po + 2r ]0 g) + V 2 / X , (80) 

19 This is shown in [34]. In [22] one can find a detailed discussion of regular solution with horizons, including a 
generalization of the Szekeres - Kruskal coordinates as well as examples of singular and degenerate horizons. 

20 An interesting exception is the integral (26) in Ref. [29] which apparently requires a more general method for its 
explanation. E.A.Davydov formulated a group theoretical approach for deriving integrals of motion in general DGS 
theories, which apparently allows to obtain more general integrals than those discussed here. 
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by solving of which we explicitly express £, rj, p, h and ip as functions of x- This is sufficient for 
finding the portrait of this physically interesting system, which will be presented elsewhere. 

To demonstrate the problems, which remain even for apparently simpler systems, we consider 
the potential U(ip) also having two additional integrals. The obvious linear integral is p = po. To 
obtain one more integral suppose that tpU^ = 2gU. Then we have the additional bilinear integral 
of the three differential equations: 

ipi] - (g + l)x + Po£ = h, X=Po + V 2 /x, rii>' = V 2 /x, ^v' = 9X- (81) 

We can exclude ip (or, rj) and thus get two equation for x an d rj (or, ip). However, this is not an 
integrable dynamical system for two functions because of its explicit dependence on £. If po = 
in the expression for the integral Iq, it can be explicitly integrated, like the previous case, but this 
gives only a 'partial' solution. This is a typical problem - bilinear integrals having ^-depending 
terms that describe a sort of a 'back-reaction' of gravity on matter. 

A detailed comparative analysis of most interesting integrable and partially integrable DGS 
systems will be given is a separate publication. In conclusion of this Section we summarize its main 
points. Dilaton gravity with scalars is in general not integrable even with formally sufficient number 
of integrals of motion. The models with massless scalars qualitatively differ from the scalaron 
models (DSG) that inevitably include non-integrability. Fortunately, in some physically important 
cases the non-integrable systems are partially integrable and therefore can be effectively studied, 
at least qualitatively. The solutions near horizons and singularities can be studied analytically, 
using exact expansions as well as iterations of the master integral equation. On the other hand, 
our approach to constructing systems with additional integrals may help to find integrable or 
partially integrable systems qualitatively close to the realistic ones. 

6 Summary and outlook 

The main new results of this paper are the following. In Section 2, the standard spherical reduction 
of the D-dimensional DVG and its Kaluza-like cylindrical reduction in the dimensions 3 and 4 are 
briefly summarized. In Section 3 we derive the transformation of the two-dimensional DVG into 
the equivalent DSG. This opens a way for applying to the vecton theory some methods developed in 
two-dimensional dilaton gravity coupled to scalar fields. In particular, the nonlinear kinetic terms 
of the vecton theory transform into completely standard potentials depending only on scalar fields 
(dilaton, scalaron, other scalars). The scalar formulation makes it easier to look for additional 
integrals of motion in the one-dimensional reductions of DSG. 21 In Section 4 we describe a more 
standard transition from the two-dimensional DVG to unified description of cosmological and static 
solutions in arbitrary gauges and write parametrization invariant Hamiltonian equations of one- 
dimensional DSG theory. The canonical constraint formalism is most convenient in searches for 
new integrals of motion, which is illustrated by the three-dimensional case. This hints at possible 
further generalizations of the approach proposed in [29] (pp. 1698-1699), where we derived two 
nontrivial one-dimensional DGS models having two additional integrals. 

Section 5 is a central part of this paper. There we introduce the simplest form of the ordinary 
differential equations describing the generic DGS theory and derive the most important properties 
of their cosmological and static solutions. In general, this theory is not integrable and we propose 
several effective approaches to its analytic and qualitative investigation. The simplest one is the 
power series expansion near the horizons and generalized expansions near the singularities (taking 
account of nonanalytic terms). We propose a more powerful approach based on an integral equation 

21 Pure dilaton gravity is a topological theory and thus reduces to the one-dimensional integrable system. There 
exist nontrivial integrable DGS models involving one massless scalar, [29]. More complex models with effectively 
massive scalar field may have one additional integral, at best, and thus generally remain non-integrable. 
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for the most important function of the system. This equation, which we call the master integral 
equation, summarizes important properties of our cosmo-stat system. In completely integrable 
cases it may be reduced to differential dynamical system. In a few instances, it was possible 
to, find a global analytic solution of the dynamical system and to draw a picture resembling the 
classical phase portraits. An example of such a portrait is presented in Appendix. 

In summary, we stress once more that a main goal of this paper is to find approaches to 
understanding the global picture of cosmological and static solutions of dilaton gravity that couples 
to nonlinear scalar fields. These theories are not integrable but, possibly, can be approximated by 
integrable ones. We have shown that there may exist additional integrals and have proposed some 
tools for exploiting this fact. Although main features of the vecton theory hint at its relevance to 
dark energy, inflation, 22 and dark matter, the precise relation is still to be uncovered, and the best 
way to solving this problem is in looking for global qualitative portraits. The difficulty is that the 
exact system is non-integrable and thus the portrait must be at least three-dimensional. 

Understanding global properties of classical solutions is also desirable for their quantization. 
The simplest approach was attempted some time ago for classically integrable gravitational sys- 
tems with minimum number of degrees of freedom (see, e.g. [44], [45], [38] and references therein). 
This primitive quantization can justifiably be criticized for not taking into account space inhomo- 
geneities, which are crucial in cosmological applications (for a review of more sophisticated ideas 
in quantum gravity see, e.g., [46], [47]). Nevertheless, we hope that such a quantization of DSG 
might be of some interest in considering properties of simplest quantum cosmological models. 



7 Appendix 

7.1 Topological portrait 

Some linear integrals of Section 5 were found in Ref.[29] for DGS with V'-independent potentials 
which, according to (70), have the integral Z(ip) tp = ryo • If also = gU, we find from (77), (78) 
one more integral (recall (68): 

ZF - g<p = p - g X = C x ■ (82) 

Then it follows that the solution can be expressed in terms of quadratures, see Eqs.(24)-(25) in 
[29]. An interesting thing not mentioned in [29] is that this integral also exists for ^-dependent 
'multiplicative' potentials U = u(ip)v(ip) when Z(tp) ip is not constant. 

Consider the second integrable model of [29], in which the potential is independent of ip and 
depends on an arbitrary function w((p). Then, in Weyl's frame V and Z are given by 



U(0 = ZV = [ 9l w(0 + g2w-\0]', 
where d£ = dip/Z(ip), = Inw(^), and evidently, 



i(0 = -(93/94) 



iTy/i+ZigM) 



[2g 3 l(0+ 9^(0} =Z, 



£/2g 3 , if 54^0. 



(83) 



(84) 



The third integral was derived by a 'brute force' in [29]. In notation of the present paper it is 

[p/ -*92h + 2g 4 % 2 \nh = C x . (85) 

The full Lagrangian and the meaning of the coupling constants can be found in the quoted ref- 
erence. As shown there, one can derive the solution in quadratures, which are not elementary if 
(74 7^ 0. If we take 54 = 52 = 0, we can obtain a simple solution expressing w as a function of h. 



22 While the relation to dark energy models is evident, the inflation was clearly demonstrated only in numerical 
studies of the massive vector meson theory in the tachyonic case m 2 < 0, [43]. 
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Introducing the scales wo and ho for w(£,) and h(£) we find the relation between normalized w 
and h, which depends only on one parameter 5 defined by the relation: 



25 + 1 = > , (86) 

where g% < for normal scalars and C\ > 0. Then the (h, u>)-portrait of our system is given by 

\h\ s 

" = IT+iiW <87) 

where e = h/\h\, —1 < h < 00, < w < 00. Now it is not very difficult to draw the picture of 
the curves describing all possible solutions. In the domain h < we have static solutions, while 
for cosmological ones h > 0. The picture looks like a phase portrait of a dynamical system in the 
(h, u>)-plane, with singular points: (0,1), (1,5), (0,0), (0, 00), (— l,oo), (00, 0). These points are 
joined by the important separating curves. The most interesting points are: the node of the initial 
singularity, (0,0), the saddle point of the horizon, (0, 1), and the most interesting 'cosmological' 
point, (1, ^), at which all cosmologies tangentially coincide. 23 

This topological portrait describes qualitative properties of static and cosmological solutions for 
different values of the parameter S, which characterizes the energy of the massless scalar field. If we 
deformed it by applying any continuous differentiable transformation preserving the singularities, 
it will represent essentially the same cosmo-static system. For example, it is easy to move all 
singularities to a finite domain of the plane. Moreover, to better understand the topological 
structure of the space of the solutions one may try to extend the portrait to the domains where, 
say, 25 + 1 <0 or h < —1, etc. We believe that the topological portrait is a most adequate global 
representation of integrable gravitational systems studied in this paper. 

7.2 Dimensions and Lagrangians 

The dimensions of geometrical and physical fields and parameters were discussed in [19]. However, 
in different papers on the affine models we used somewhat different conventions and notation. Here 
we first give several comments on these matters. In the pure geometrical part, there is no problem 
at all as we have only the dimension of length. Thus [sij] = [ajj] = L~ 2 , [I] = L°, geometric 
Lagrangian (1) is of dimension L~ D , and geometric action is therefore dimensionless. In transition 
to a physical picture, we define dimensionless tensor densities g y , P J by varying Lagrangian (1). 
To make the densities dimensionless we multiply it by a constant 27 where [7] = L D_2 : 

d daij 

Using this definition it is not difficult to prove the following identity (/ is a dimensionless constant): 
|A fl | = I det(g« + / F)| = 7 D I det( Sij + T 1 ay)| {D ~ 2)/2 • (89) 

With this identity, one can find the conjugate Lagrangian density C* = £*(g u , f 3 ): 

F)C* r)C* 
= ^ (" Vf) " 5 »r, —j. a l3 = —j . (90) 

This result was apparently known to Einstein in the special (four-dimensional) case, when v = 1/2 
and I = 1. Following his approach we find the 'physical' Lagrangian (2) as described in [18] - [22]. 



23 Even more interesting but much more complex is the portrait of the model with massless vector (when 32 ^ 0) 
and scalar, see Eq.(30) in [29]. We will discuss it in a separate publication. 
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Eqs.(88) - (90) play the key role in this derivation. When considering the physical Lagrangian 
it is better to use the gravitational constant k = G/c A to define the physical fields A4 and F^ 
proportional to Oj and fij but having the standard dimensions and satisfying the relation [kA 2 ] = 1, 
[nFl] = \T 2 , [19]. In this paper we choose the units c = k = 1 and thus A{ is dimensionless and 

is of dimension L _1 . In fact, we use notation a. L for the dimensionless potential in this system 
of units and fij for the corresponding field tensor = dij — a,jj ; this means that [A] = L and 
[A 2 A] = 1. In Eq.(2) and in the body of this paper we always apply this agreement on notation. 
As we do not use the tensor densities in the main text, the bold-face notation is used only to 
fij fi = f 2 and a, a* = a 2 . 

One final remark on the dimensions of constants and fields. In our units, we have essentially 
three dimensional constants of geometric origin: A, A, fx = m 2 (k u and 7 are auxiliary dimensional 
constants having no relation to affine geometry). Coordinates t,r,u,v are of dimension L, the 
main fields a^, if), 99 are dimensionless while [q] = L t 1 . When we denote the scalaron by ip (in 
Section 5) we silently make it dimensionless by including m in its definition, ip = q/m. 

7.3 Three-dimensional pure vecton theory 

Let us consider the vecton Lagrangian (42) in the LC gauge (when = 7 and 2a = ln\h\ = F) 
and without transition to the scalaron description 

— fj,(pa e — 2A<ph, (91) 

where h = ee~ F and Ao = A 2 A are dimensionless while \x = m 2 is of dimension L~ 2 . The 
Hamiltonian constraint that must be derived before passing to LC coordinates is 

•H« = X Q{p h 2 /h - ipF + ^ipal + 2kiph = (92) 

while other equations are simply derived by varying Lagrangian (91). Instead of the Hamiltonian 
equations we write the equations of motion in the first order form similar to (70). Defining 
£ = (ip 2 /2), (ipa)/h = b,A = a/%, \i\ = fi/Xo, and denoting d/d^ by the prime we find: 

b' = -n x A , H' = fiA 2 H , 2 fa' = bH , & = (2Af + A b 2 /2) H , (93) 

where we omitted index e of a and A. In the massless limit fx = [i\ = and therefore H = Hq, 
b = 60 are independent of f . This immediately gives the three-dimensional Maxwell - Einstein 
solution. Indeed, we see that 

HO = H X (0 , a = a + (b H /2) In , X = Xo + (A b 2 H /2) In + 2A(f - fo) , (94) 

and thus returning to the variable ip we find (taking Xo = and denoting Ao 60-^0 = c): 

h = H £ dt(j^ + 2A ] j = H 1" dp ^ + 2A^ . (95) 

Equations (94)- (95) define the general solution of the pure dilaton gravity with the potential 
proportional to the expression in the brackets, see Eq.(67). 

To find approximate solutions for small but finite values of fx, \x\ we may take H = Hq, b = bo, 
expressions (94) for cto(f), Xo(£)j ^0 = «o/xo and thus construct successive approximations for 
bi, Hi,cii,Xi, Ai = (Xi/xi- Our equations (93) or their integral formulation are convenient for 
deriving the solution near horizons and near singularities at f — > 0,oo. In fact, they are simpler 
than the approximate equation studied in [18] - [21]. At the horizon at f = fo> when ao = Xo = 
and .ffo,^4o are finite, we have the convergent expansions of the solution in powers of f — fo- 
Moreover, one can see that the radius of convergence cannot be larger than fo (it can be lesser 
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if there is a second horizon). Near singularities we can also find expansions in powers of £ if we 
correctly take into account the logarithmic terms seen in the zeroth approximation (95). 

A better though much more complex iterations can be derived with the aid of the integral 
equations for A(£) similar to (73), (74), if we take the zeroth approximation A = Aq, solve the 
equations for H and b, solve the resulting equation for A(£), and then repeat this cycle. A simpler 
approach is to first find and use an additional approximate integral. On this way, the scalaron 
formulation looks more promising. For example, we have seen that the potential for D = 3 
depends on q 2 /^ and there exists the additional integral (64), which presumably can be used for 
constructing more effective approximations. 
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